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Abstract 

We use QCD sum rules to determine the difference between moments of the non- 
singlet structure functions. This combination decouples from the singular behaviour of 
the structure functions near x = 1 as calculated in the quark-gluon basis and thus should 
lead to improved sum rule predictions. However, we find there are still very large errors 
due to higher order corrections. In order to refine the error analysis, we study the effect 
of renormalon ambiguities on our QCD sum rules results. 
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1 Introduction 

The application of QCD sum rules to deep inelastic scattering has been pioneered by Belyaev 
^ ; and Ioffe [|l], || and subsequently developed further by Balitsky et al |J. In 0] we used QCD 
sum rules to determine the leading twist-two non-singlet operator matrix elements (OMEs) 
^h. contributing to deep inelastic scattering from nucleon targets and an error analysis was carefully 
performed. The comparison with the measured values was disappointing and we investigated the 
source of errors in detail. We found large sources of error from both the neglect of higher order 
terms in the expansion in powers of 1/p 2 where p is the momentum flowing through the nucleon 
source ( "higher-dimension" terms) and the neglect of higher order terms in the perturbative 
expansion in powers of the QCD coupling. The phenomenological analysis is very sensitive 
to the singular contributions in the structure functions at x — 1 corresponding to unphysical 
quark-gluon kinematics. The authors of || have suggested that this problem can be evaded if 
one uses the so-called "Ioffe-time coordinate space distributions", as a suitable alternative to 
the conventional parton momentum distributions. Here we examine another way of suppressing 
the singular contributions at x — 1. To do this, we consider specific combinations of moments 
which eliminate the effects of the singular corrections at x — 1. Unfortunately, there is not a 
great improvement in the QCD sum rule predictions. Again there seem to be anomalously large 
corrections at O (u s (qq) 2 ' j \p 2 ) 2 ) coming from the radiative correction to the quark condensate 
term. We also investigate further corrections due to renormalon ambiguities giving power law 
corrections which cannot be absorbed in the terms arising from quark and gluon condensates. 
Indeed, renormalon singularities and their implications for theoretical predictions have been 
studied for some time and it has now become clear that one cannot sensibly talk of power- 
suppressed corrections if the ambiguities in the leading term are not under control [§, 0, §]. 
The reason the renormalon ambiguity arises is that the QCD sum rule provides an estimate for 
an OME where the operator is normalised at the scale p 2 . Estimating the OMEs' dependence 
on p 2 involves perturbative calculations that introduce the renormalon ambiguities. There is a 
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simple correspondence between renormalon positions and the power corrections to fixed-order 



perturbative predictions evaluated with an infrared cutoff [10|. Following Webber's approach 
|TH , we impose such a cutoff by introducing a small mass m 2 into the gluon propagator and 
examine the effect on the QCD sum rules determination of the structure functions moments and 
whether it is process dependent or not, i.e. whether the effect is important for some moments 
while it is negligible for others. 



2 QCD sum rules for the difference of moments 

Our starting point is the result obtained in |Q for the quark distributions (valid for intermediate 
values of x) : 
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Eq(z) = 1 — e~ z , E2{z) = 1 — e~ 2 (l + z + ), W is the continuum threshold while L = 
ln(M/A)/ln(/i/A) takes account of the anomalous dimension of the currents. The dependence 
of Afy on the Borel parameter M 2 is expressed via the mass sum rule ||12|]p| 
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In equations (H) and (||]), the first term on the LHS comes from the resonant saturation with 
a nucleon double pole while the second parametrises some of the effects of the non-resonant 



^dn the sign of the anomalous dimension of baryonic currents is wrong resulting in a wrong expression 
for \ 2 N {M 2 ) in Q]. Taking account of this fact affects only slightly the results and conclusions of [Q. 
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background corresponding to a single pole term. The RHS comes from an evaluation of the 
deep-inelastic scattering process in perturbative QCD. 

In ||, these relations were used to estimate the quark distribution moments M% defined by 

Ml = f 1 dxx n u v (x,Q 2 ) 
Jo 

M? = [ dxx n d v (x,Q 2 ) (4) 
Jo 

which are related to the OMEs A^ ,d = (iV|0^' d | N) and the Wilson coefficients C n in the operator 
product expansion (OPE) via 

M n (Q 2 ,a s (Q 2 )) = J2AM,a s (i4))C(j\Q^ f il,a s (i4)) (5) 

T 

where the sum runs over the twist r of the operators O n and fi^ is the scale at which the 
operator O n is renormalised. 

A phenomenological study was performed for the OMEs and it revealed a large discrep- 
ancy when comparing with experiment. In j|| we argued that this discrepancy was due to a 
breakdown in the perturbative analysis used in the QCD sum rules. In particular, the O (a s a 2 ) 
radiative corrections give contributions to the structure functions proportional to tjz^) (which 
must be regulated at x = 1) and these terms gave very large contributions to the moments 
suggesting that the perturbation series in a s does not converge. Further, there are terms pro- 
portional to 5(1 — x) (not displayed in eqs ([]]) and (|2[) which apply at intermediate x) which 
give very large higher dimension terms in the moments and suggest the expansion in powers of 
1/p 2 does not converge either. 

One point of view promoted by Belyaev and Ioffe is that the QCD sum rule predictions 
should be used only for the structure functions at intermediate x, keeping away from the 
troublesome singularities at x = 1. However we find this unconvincing as such terms come 
from unphysical quark and gluon singularities which must be averaged over when making 
comparison with data. This averaging procedure will feed the effects of the singular terms 
at x = 1 to intermediate x causing the same discrepancy with data as is found in the moment 
relations (which represent a particular form of averaging). Another possibility that has recently 
been proposed || is that the problem can be evaded if one uses the so-called "Ioffe-time 
coordinate space distributions" , as a suitable alternative to the conventional parton momentum 
distributions. 

In this paper, we wish to explore a more straightforward way of eliminating the singularities 
at x = 1 for it is easy to do so even in the conventional picture. In order to achieve this, one 
considers the difference of two consecutive moments. To see that this eliminates the effect of 
singularities at x = 1 we note that the structure function has the form (c.f. eqs (||) and (||)): 

q{x) = h{x) + f 2 (x)6{l-x) + ^ (6) 

1 — x 

where fx, f'2 and f'3 are non-singular functions. By integrating, one obtains the moments: 

M n = f 1 x n q(x)dx = C f x {x)x n dx + / 2 (1) + f 1 ^-x n dx. (7) 
Jo Jo Jo 1 — X 

As mentioned in ||], in order to find the "physical" moments one needs to regulate the third 
term: 

A n = f 1 h(x)x n dx + / 2 (1) + f 1 h{x) f ~ dx . (8) 
Jo Jo 1 — X 
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(a) (b) 

Figure 1: Twist-two non-singlet OMEs' difference following from QCD sum rules together with 
the error estimates. Also shown are the experimental determinations of the OMEs' difference 
(diamonds). 

Now we see that if we subtract A n from A n+ i the effect of the singularity 5(1 — x) disappears 
while that of is reduced to an integral of a regular function (the effect of the regulating 

term cancels). 
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Using eqs (HD and (0) we may write down the QCD sum rules for this difference in moments 
and arrive at: 
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(a) (b) 

Figure 2: Analysis of the importance of the various contributions to the QCD sum rules to 
the OMEs' difference. In (a) we consider the up quark operators. The stars show the results 
neglecting the gluon condensate and the 0(a s a 2 ) corrections. The dots show the results includ- 
ing the gluon condensate. The triangles show the complete result, (b) shows the equivalent 
quantities for the down quark. 
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The same phenomenological analysis as in 0] can be done for these sum rules and the results, 
for the same numerical values of the parameters as in JI], (a = 0.55 ± 0.20 GeV 3 , b = 0.45 ± 
0.10 Gel/ 4 , fi = 0.5 GeV, /i 2 , = 1 GeV 2 , A = 125 ± 25 MeV, m = 1.00 ± 0.15 GeV and = 
2.3 ± 0.23 GeV 2 ) are shown in Figs. and @. 

One may see from Fig. [I] that there is still a discrepancy between the prediction of QCD 
sum rules and the experimental measurement. To analyse why, we show in Fig. |^ the effect of 
including various contributions. One may see that the O (ct s a 2 ) contribution is still very large 
suggesting that the perturbative series is not convergent. 
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Figure 3: Diagrams involved in the calculation of the OME. The solid lines represent quarks, 
the wavy lines represent gluons and the circle represents the operator. 

3 Renormalon effects 

Given that the QCD sum rule predictions fail to reproduce well either the moments or the differ- 
ence of moments it is perhaps appropriate to consider other possible sources for the discrepancy 
beyond the immediate (and depressing) possibility that the expansion in higher dimension terms 
and higher order radiative terms is not convergent. We have identified one possible further error 
coming from renormalons. 

To see how renormalons affect our analysis of the moments, we look again at the OPE form 
for the moments, eq @. There are two scales in our QCD sum rules analysis: p 2 and Q 2 . The 
sensitivity to nucleon momentum scale p 2 , which we assume to be large and euclidean in order 
to justify the perturbative evaluations of the graphs, resides in the OME A n in eq fl5|), while 
the dependence on the momentum transfer scale, Q 2 , resides in the coefficient function C n . 
By going to large Q 2 the coefficient functions may be reliably calculated, as discussed in [|]]. 
However, the p 2 dependence cannot be reliably calculated in perturbation theory because the 
matching condition of QCD sum rules requires values of p 2 K,l GeV 2 where non-perturbative 
corrections are important. Once the Q 2 -dependence is removed using the OPE the only scale 
other than p 2 is fi 2 , the operator renormalisation scale. The latter is in principle arbitrary, 
the operator /^-dependence being cancelled by the /^-dependence of the coefficient functions. 
In practice, however, in evaluating the OME by sum rules, the normalisation scale must be 
chosen of the same order as the only other scale p 2 to avoid the appearance of large radiative 
corrections (which have not been calculated and thus do not show up explicitly in the analysis). 
Thus, one should interpret the result of the QCD sum rule predictions as being for the operator 
renormalised at a scale p 2 . Then, in order to calculate the OME at a fixed scale fi 2 appropriate 
to the comparison with the matrix elements extracted from deep inelastic scattering processes, 
one must determine the /^-dependence of the OME A n (n 2 ) 



(N(p)\d Pai ... an (fi 2 )\N(p)) = A n ^ 2 )pPp a \..p an +••• • (12) 

to allow A n (p 2 ) to be related to A n (fil). 

To estimate the /^-dependence of the OME, one should evaluate the graphs in Fig. [j| 
The "renormalon" contributions are described by graphs such as those in Fig. [| and in order 
to take account of their effects, one should introduce an infrared cutoff in the integrals to be 
done when evaluating the graphs. Equivalently one may introduce a mass m g into the gluon 

propagator and expand in powers of 37%, a procedure which has the advantage of maintaining 



Lorentz invariance |TT[ • Using this technique after some calculation one arrives at the following 
results: 
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Figure 4: Diagrams leading to renormalon ambiguities in A n . The solid lines represent quarks 
and the wavy lines represent gluons. The circles inserted in the gluon lines refer to fermion or 
gluon loops and the sum should be over all such insertions. 
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We now introduce the renormalisation scale fi and we get, up to order g , 
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We see from this equation that the relation between matrix elements for operators renormalised 
at different scales involves corrections of O corresponding to the mixing of operators of 



different twist. This difference gives a further source of error in the QCD sum rule predictions 
which we now consider. As we argued above, one should interpret the QCD sum rule predictions 
as giving A n {p 2 ) and so the left hand side of the equation has the form ^"^2 for the nucleon 

pole contribution. The ln(^) and the ln(q^) terms give, when Borel transformed, terms 
proportional to ln(^^) and ln(^) respectively. These terms can be neglected if, for the 

purpose of getting a rough estimate of the renormalon effects, we choose m 2 ~ | A*o I — ^GeV 2 to 
be in the QCD sum rules' "overlapping window" (M 2 ~lGeV 2 ). For this choice (m 2 ~ - 1,1 5 

2 2 2 

can also ignore the In 3^ term. That leaves us with the terms proportional to + -jpr ■ The 
former contributes additional terms corresponding to different poles on the phenomenological 
side of the sum rule 

1 =1 =± A. 

m 4 I m? 1 m4 /-1 oA 
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—p 2 (p 2 — m 2 ) 2 p 2 (p 2 — m 2 ) 2 p 2 — m 

The last term corresponds to a single pole contribution so it can be absorbed by the unknown 
non-resonant background effect. The first term is also expected to be small since the structure 
of eq ( |TTD applies only for p 2 Euclidean and far from 0. To enforce this we introduce a threshold 
cutoff in the Borel integration of this term, leading to a vanishing contribution. Finally, the 

2 

contribution of the second term is of the same form as that of the term proportional to —§■ in 
eq (fT7P. Combining the two with — /1q~ + m 2 ~lGeV 2 gives 

A n (p 2 = 0(ti)) * ^)(l + 2^F(n)^). (19) 

Physically, we view the gluon mass m 2 as an infrared matching parameter which represents 
the scale below which we switch from the perturbative to the non-perturbative domain and 
which is in principle much greater than Aqcd- We shall take it to be a non-perturbative 
and process-independent parameter (c.f. |TJJ). If we denote by A® CDSR the QCD sum rules 
predictions of M then, to first order in as, we have 



as{l4) C F 
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2^ 



An(d) = Ar DSR (1 - 2^g^ F{n) 'jL j . (20) 

4 Discussion 

Eq (|20|) shows that renormalon effects lead to a correction of the QCD sum rule estimate 
of the OME as measured in deep inelastic scattering. Following fllTf , we shall assume that 



the n-dependence given by the perturbative calculation with a gluon mass correctly gives the 



mi 



relative magnitude of the renormalon corrections, the only unknown being the ratio —§-. We 

Mo 

compare this renormalon effect with the deviation of the A® CDSR from the experimental 
data obtained by fitting the relation 

A**>(fi) = A% CDSR (l + di) (21) 



between our QCD sum rules predictions A® CDSR and our experimental data A'^ lta {ji^) for the 
OME. The results are shown in Table |1] for the u-quark and Table for the d-quark where we 
have taken the values " s ^°~ 2 ~ 1 GeV ^ ~0.054, — 1 GeV 2 and m 2 ~l GeV 2 . 
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n 


2 


3 


4 


5 


6 


7 


8 


9 


F(n) 


0.5 


-0.83 


-2.75 


-5.12 


-7.85 


-10.89 


-14.20 


-17.75 


j^data (r\\ 


0.11 


0.043 


0.021 


0.011 


0.0064 


0.0039 


0.0025 


0.0017 




0.32 


0.31 


0.32 


0.34 


0.35 


0.36 


0.37 


0.38 


d A (n) 


-0.67 


-0.86 


-0.94 


-0.97 


-0.98 


-0.99 


-0.99 


-0.995 


Ren A (n) 


0.072 


-0.12 


-0.40 


-0.74 


-1.13 


-1.57 


-2.04 


-2.56 


d A (n)/Ren A (n) 


-9.26 


7.18 


2.36 


1.31 


0.87 


0.63 


0.49 


0.39 




-0.062 


-0.023 


-0.0098 


-0.0047 


-0.0025 


-0.0014 


-0.0008 






-0.0035 


0.0088 


0.012 


0.013 


0.013 


0.012 


0.011 




d K (n) 


16.49 


-3.54 


-1.78 


-1.35 


-1.19 


-1.11 


-1.07 




Ren K (n) 


17.04 


-10.20 


-9.58 


-11.22 


-13.58 


-16.34 


-19.46 




d K (n)/Ren K (n) 


0.97 


0.35 


0.19 


0.12 


0.088 


0.068 


0.055 





Table 1: Results for the u-quark for the numerical values — IGeV 2 and m 2 g ~lGeV 2 . 
K(n) = A(n + 1) - A(n); A data (n) and K data {n) are the experimental data while A QCDSR (n) 
and K QCDSR (n) represent the QCD sum rules predictions. d A (n) represents the deviation 
between data and QCD sum rules predictions, Ren A (n)= —2 as ^ C F F(n)^§- represents the 

"renormalon" correction, and d A (ii)/Ren A (ji) represents the proportion of these two quantities. 
d K (n), Ren K (n) and d K (n) /Ren K (n) represent the same quantities but for the difference of 
moments K(n). 



n 


2 


3 


4 


5 


6 


7 


8 


9 


j^data (yy^ 


0.031 


0.011 


0.0047 


0.0023 


0.0012 


0.0007 


0.0004 


0.0003 


A QCDSRf n \ 


0.70 


0.61 


0.55 


0.50 


0.47 


0.44 


0.41 


0.39 


d A (n) 


-0.95 


-0.98 


-0.991 


-0.995 


-0.997 


-0.998 


-0.999 


-0.999 


Ren A (n) 


0.072 


-0.12 


-0.40 


-0.74 


-1.13 


-1.57 


-2.04 


-2.56 


d A (n)/Ren A (n) 


-13.25 


8.18 


2.50 


1.35 


0.88 


0.64 


0.49 


0.39 




-0.020 


-0.0064 


-0.0024 


-0.0011 


-0.0005 


-0.0003 


-0.0001 






-0.087 


-0.061 


-0.045 


-0.036 


-0.029 


-0.025 


-0.021 




d K (n) 


-0.77 


-0.89 


-0.95 


-0.97 


-0.98 


-0.99 


-0.995 




Ren K (n) 


1.42 


2.37 


3.40 


4.43 


5.44 


6.43 


7.38 




d K (n)/Ren K (n) 


-0.54 


-0.38 


-0.28 


-0.22 


-0.18 


-0.15 


-0.13 





Table 2: The same explanation as for the previous table but for the d-quark. 
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One can do the same analysis for the difference of moments K n = A n+ i — A n , to obtain 
K ,»s jv-qcdsr (, _ ^l)C F A% DSR F(n + 1) - A*°DSR F{n) 

^nKH-OJ — \ ^7T j^QCDSR _ j^QCDSR ^2 J 

Kfiti) = K% CDSR {l + d«) (23) 

where K® CD denotes the QCD sum rules predictions for the difference of moments K n , and 
j^data are correS p 0ri di n g experimental values. These results are also shown in Tables |TJ 
and 0. 

We see from both these tables that, for higher moments (n>5), the magnitude of the cutoff- 
dependent "renormalon" contribution for the OME A n is comparable to, or greater than, the 
experimental deviation d^. This situation extends even to lower moments in the case of the 
difference of moments K n . This means that renormalon corrections cannot be neglected in 
the QCD sum rule evaluation of the higher moments. It is possible to account for some of 
the discrepancy between the QCD sum rules and experiment via the renormalon term through 
a suitable choice of m 2 . However, at best this only alleviates the problem slightly. A more 
cynical view is that the new corrections identified here are just another source of large higher 
dimension corrections to QCD sum rules which render the whole method of doubtful use in 
determining the nucleon properties. As we have seen, the problem persists for the subtracted 
moments so we cannot blame the unphysical quark and gluon singularities at x — 1 for the 
failure of the sum rules. The disappointing conclusion of this is that in determining the leading 
twist OMEs, where good experimental measurements are available to check the results, the 
QCD sum rule method fails largely because of uncontrollable higher dimension corrections. In 
our opinion, this casts doubt on the QCD sum rule predictions for other nucleon properties 
such as higher twist matrix elements (where the troublesome higher order corrections have not 
been computed). 

Acknowledgement: We are very grateful to M. Birse for valuable discussions. 
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